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ENUMERATIVE COMBINATORICS OF SIMPLICIAL AND CELL COMPLEXES: 
KIRCHHOFF AND TRENT TYPE THEOREMS 


SYLVAIN E. CAPPELL AND EDWARD Y. MILLER 

ABSTRACT 

This paper considers three separate matrices associated to graphs and (each dimension of) cell complexes. It relates 
all the coefficients of their respective characteristic polynomials to the geometric and combinatorial enumeration of 
three kinds of subobjects. The matrices are: the mesh matrix for integral d-cycles, the mesh matrix for integral 
d-boundaries, and the Kirchhoff matrix, i.e., the combinatorial Laplacian, for integral (d-l)-chains. 

Trent’s theorem states that the determinant of the mesh matrix on 1-cycles of a connected graph is equal to the 
number of spanning trees |29l |30). Here this theorem is extended to the mesh matrix on d cycles in an arbitrary 
cell complex and, new even for graphs, to enumerative combinatorial interpretation of all of the coefficients of its 
characteristic polynomial. This last is well defined once a basis for the integral d-cycles is chosen. Additionally, a 
parallel result for the mesh matrix for integral d-boundaries is proved. 

Kirchhoff’s theorem states that the product of the non-zero eigenvalues of the Kirchhoff matrix, i.e., combinatorial 
Laplacian, for connected graphs equals n times the number of spanning trees with n the number of vertices. Lyons has 
generalized Kirchhoff’s result on the product of the non-zero eigenvalues of the Kirchhoff or combinatorial Laplacian 
on (d-l)-chains to cell complexes for d > 1 [18| . The present analysis extends this to all coefficients of the characteristic 
polynomial. 

An evaluation of the Reidemeister-Franz torsion of the cell complex with respect to its integral basis gives relations 
between these combinatorial invariants. 
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1. Introduction 


This paper considers three separate matrices associated to graphs and (each dimension 
of) cell complexes. It relates all the coefficients of their respective characteristic polynomials 
to the geometric and combinatorial enumeration of three kinds of sub-objects. Their de¬ 
terminants and more generally their characteristic polynomials are interpreted in terms of 
correspondingly three different aspects of the enumerative combinatorics of the graph or cell 
complex (in each dimension). 

The oldest such matrix, of the three studied here, is Kirchhoff’s matrix, which is the com¬ 
binatorial Laplacian A 0 unim on 0-chains. For graphs, A 0 is the degree matrix minus 
the adjacency matrix. For a connected graph with n vertices, Kirchhoff showed that the 
product of the non-vanishing eigenvalues of A 0 is equal to n times the number of spanning 
trees. This interpretation of the product of the non-vanishing eigenvalues of the combina¬ 
torial Laplacian on (d-l)-chains was generalized to higher dimensional complexes by Lyons 
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[18j as an explicit weighted sum over pairs of combinatorial objects which he called bases 
and cobases, and here (consistent with the classical terminology for graphs) are called (d- 
dimensional) spanning forests and ((d-l)-dimensional) spanning coforests. In the restricted 
context of d-dimensional complexes which are rationally acyclic in dimensions less than d, 
parallel theorems to those of Lyons on the product of the eigenvalues of the Krichhoff matrix 
had earlier appeared in the work of Duval, Klivans, and Martin |12j . 

In the case of graphs the complete characteristic polynomial of the Kirchhoff matrix had 
been studied and each of its coefficients had been evaluated by an enumerative combinatorial 
formula, a sum over special subgraphs, forests with specified number of edges with associated 
weights For cell complexes, analysis of the characteristic polynomials is developed here 

giving an analogous interpretation for all coefficients of the polynomial in each dimension. 

The second type of matrix studied here is the higher dimensional analogue of the mesh 
matrix introduced by Trent pifJl [30] to study enumerative combinatorics of graphs. It is 
based on the idea that the integral cycles have a natural pairing to the integers, yielding 
a mesh matrix once a basis is chosen. As shown by Trent, the determinant of the mesh 
matrix of a connected graph equals the number of spanning trees. Here this classical work 
is extended in two new directions. One extends the definition and results about the Trent 
matrix for graphs to simplicial or cell complexes, obtaining a mesh matrix in each dimension 
d with determinant a sum over (d-dimensional) spanning forests with explicit weights. Next, 
we give a geometrical interpretation for each coefficient in the characteristic polynomial of 
the mesh matrix of graphs as well as in the higher dimensional context. The mesh matrix 
is defined for each choice of integral basis for the d-cycles; its determinant is independent of 
this basis, while the coefficients of other powers of t in the characteristic polynomial depend 
on this choice. We see that in §5, the coefficient of each power of t is shown to be sum over 
geometrical objects, (d-dimensional) k augmented spanning forests, with positive integer 
weights depending in a specified way on the chosen basis. This is of interest and novel even 
in the graph setting. Moreover, these weights are here given for cell complexes and graphs 
several topological and combinatorial interpretations in settings of geometrical interest. 

Thirdly, the integral boundaries of chains in each dimension also have a natural pairing 
to the integers. Once an integral basis is chosen a matrix arises. Its determinant is the 
sum over (d-dimensional) spanning coforests with explicit positive integral weights which 
are studied below. Again each coefficient of powers of t in the characteristic polynomial 
has an expression as a sum uniformly over (d-dimensional) k reduced coforests with certain 
positive integral weights. 

The geometrical objects being enumerated with weights in the formulae for each of the 
three settings are thus different. Let rid = dim Cd(X; R), Zd = dim Z d (X;R), b d = 
dim Bd(X;R), the dimensions of the d-chains, d-cycles, d-boundaries. For the reader’s 
convenience, these are tabulated in the following table. 


ENUMERATIVE COMBINATORICS OF SIMPLICIAL AND CELL COMPLEXES: KIRCHHOFF AND TRENT TYPE THEOREMS 



for connected graphs (d=l) 

for cell complexes 

Matrix 

Types 

product 

non-vanishing 

eigenvalues 

coefficient of t L 
of characteristic 
polynomial 

product 

non-vanishing 

eigenvalues 

coefficient of t L 
of characteristic 
polynomial 


Geometrical objects summed over to yield algebraic invariant 

Kirchhoff’s 
Combinatorial 
Laplacian 
on (d-l)-chains 
rid -1 x rid -i 

spanning 

trees 

L= 1 

forests of size m 
L = no — m 

pairs, spanning 
forest and 
spanning 
coforest in it 

L = rid -i — bd -i 

pairs, forest 
of size m and 
spanning 
coforest in it 

L = n d ~ i - m 

0 < m < b d -i 

Trent 

Mesh matrix 
on d-cycles 

Z d X Zd 

spanning 

trees 

L = 0 

k augmented 
spanning forests 
L = k 

spanning forests 

L = 0 

k augmented 
spanning forests 

L = k 

0 < k < z d 

Mesh matrix 
on d-boundaries 

bd X bd 



spanning coforests 

L = 0 

k reduced 
spanning coforests 

L = k 

0 < k < bd 


In the Kirchhoff case the pairs, (forest, spanning coforest in it) are in dimensions (d,d-l) respectively; in the cycle 
mesh matrix case the k augmented spanning forest is in dimension d; in the boundary mesh matrix case the k reduced 
spanning coforest is in dimension d. A k augmented spanning forest is a subcomplex for which removing some choice 
of k d-dimensional cells yields a spanning forest; a k reduced spanning coforest is a subcomplex for which adding 
some choice of k d-dimensional cells yields a spanning coforest. 

Additionally, for the d-cycles given a spanning forest there is a geometrically defined basis 
for the rational cycle group Z d (X\ Q)\ again for the d-boundaries given a spanning forest in 
dimension d + 1 there is a rationally defined basis for the rational boundary group B d (X] Q ). 
Consequently, one obtains two types of geometrically defined mesh type matrices. In section 
7, in each of these cases the complete characteristic polynomial is evaluated combinatorially. 

A spanning tree for a connected graph with n vertices is a choice of n — 1 edges which 
together carry no 1-cycles. More generally, even when the graph is not necessarily connected, 
a choice of spanning tree for each connected component is called a spanning forest. In this 
possibly disconnected context, Trent’s theorem for graphs asserts that the determinant of 
the mesh matrix for cycles is the number of spanning forests [29, [30]. Correspondingly, 
the generalization to the complete characteristic polynomial below is expressed in terms of 
spanning forests and its extension, k-augmented spanning forests. 

The mesh matrix of 1 cycles for a graph and more generally for d cycles of a CW complex 
X is easily explained. The d chains, C d (X; Z ) have a basis of oriented d cells, and this yields 
a combinatorially defined pairing by setting < crj,ak >= 0 if the cells are different and +1 
if j = k. By linearity this gives a pairing of integral chains 

<•,■>: C d (X;Z)xC d (X;Z)^Z. 

By definition the integral d cycles, Z d (X;Z) are the kernel of the boundary mapping d : 
C d (Z]Z) —> C d _i(X\ Z). Now including the d cycles Z d (X]Z) C C d (X;Z) the composite 
Z d (X- Z) x Z d {X] Z) C C d {X] Z) x C d (X ; Z) —> Z gives the pairing <.,.>. 

Taking a basis for the d cycles, which is a free abelian group, say {z r }, there is the 
associated mesh matrix on d cycles 

M{z r } = {< Zj,z k >}• 
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An analogous matrix is associated to the integral d boundaries, B d (X] Z) = image of d : 
Cd+i(X] Z ) — > C d (X ; Z ) giving a mesh matrix for d boundaries, once a basis is chosen for 
the integral boundaries. 

In this paper, the determinant of the mesh matrix on d-cycles is shown to be a simple 
positive integer weighted sum over what are called here “spanning forests” in a finite cell 
complex A", extending the classical notion of spanning forests for graphs. These spanning 
forests are subsets of d cells maximal with respect to the property of carrying no d cycles. 
This extends Trents’s theorem to each dimension, say d, of a cell complex. Two distinct 
combinatorial evaluations of the weights are here provided. 

There is a corresponding notion of k augmented spanning forest which are just spanning 
forests with k added d cells. Having picked a basis for the integral d-cycles, the mesh matrix 
on d-cycles is a well-defined integral valued z d x z d square matrix. It is shown that the 
(z d — k) th elementary symmetric function of the eigenvalues of of the mesh matrix for cycles 
is a sum over corresponding k augmented spanning forest with positive integral weights. 
These weights are given an explicit dependence on the basis of d-cycles chosen. This extends 
Trent’s theorem to the complete characteristic polynomial of the mesh matrix for cycles. 
The weights are given in two quite different formats (but equal). This is of interest even in 
the graph, d — 1, case. These definitions and results appear in (J2J ( The spanning forests 
are called bases in Lyons [TH] . ) 

Additionally, the mesh matrix for d boundaries has a parallel development for cell com¬ 
plexes; although it is trivial for graphs. Having chosen a basis for the integral boundaries, it 
is a well defined matrix. The determinant is obtained by a positively integer weighted sum 
over spanning coforests. Having picked a basis for the integral d-boundaries, the (b d — k) th 
elementary symmetric function of the eigenvalues of the mesh matrix for d boundaries is 
shown to be a sum over k reduced spanning coforests with appropriate positive integral 
weights depending on the choice of basis. A k reduced spanning coforest is just a spanning 
coforest with k cells removed. These definitions and results appear in [J5] Again two distinct 
evaluations of the weights are provided. 

In section 7, for each of two mesh matrix cases, the complete characteristic polynomial is 
evaluated combinatorially in yet another way, using a geometrically defined basis of the ra¬ 
tional d cycles and boundaries. This employs a fixed choice of a spanning forest in dimension 
d and d + 1 respectively. 

In the Kirchhoff case and in the mesh cycle and mesh boundary case there is a formal 
analogy. The combinatorial Laplacian on d chains, the mesh matrix for d-cycles, and the 
mesh matrix for d boundaries each are represented by matrices of the form, A ■ A*. The 
Cauchy-Binet theorem applies as an algebraic step in the analysis of such A ■ A 1 . The non¬ 
zero eigenvalues of A- A t are real, non-negative numbers, their positive square roots are often 
referred to as the singular values of the matrix A. 

This algebraic procedure applied in the combinatorial Laplacian setting recapitulates the 
approach of Lyons who proved that the product of the non-zero eigenvalues of this Laplacian, 
acting on (d-1) chains, equals the a sum over (d dimensional) spanning forests and ( (d-1) 
dimensional) spanning coforests contained in them with positive integral weights in the 
terminology used here [bases and cobases in his treatment] [18] . Here this result is extended 
to the coefficients of the complete characteristic polynomial using a sum over pairs, a (d 
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dimensional) forest and a ((d — 1) dimensional) spanning coforest in it. These results appear 
in §8. 

There is a large and growing literature on the general topic of generalizations, analogues 
and applications of Kirchhoff type theorems. A small sample of relevant combinatorially 
oriented papers and monographs would include (1847) G. Kirchhoff [17], (1889) A. Cayley 
[7], (1935) H. Whitney [33], (1954) H. M. Trent [29], (1955) H. M. Trent [30], (I960) N. 
I . Biggs [2], (1961) A. Nerode and H. Shank [2i] , (1967) P. R. Bryant [4], (1971) W. K. 
Chen [ID], (1971) H. Shank [28], (1976) S. B. Maurer [ID], (1976) J.W. Moon [22], (1978) 
S. Chaiken and D. Kleitman [8], (1979) W. Tutte [31], (1982) S. A. Chaiken [9], (1983) G. 
Kalai [2], (1992) R. M. Adin [T], (2000) R. Kenyon [T5], (2009) A. Duval, C. Klivans, and 
J. Martin [12] . (2009) R. Lyons [TS], (2009) A. Petersson, (2011) R. Kenyon [TR], (2013) M. 
Catanzaro, V. Chernyak, and J. Klein [6j, the items after 1983 in higher dimensional settings. 

The determinants of the mesh matrices for d cycles divided by that for d boundaries is 
directly related to the determinant of a corresponding matrix for the homology. A compu¬ 
tation of the Reidemeister-Franz torsion of the CW complex A", yields interesting relations 
with the above combinatorial items. These definitions and combinatorial identities appear 
in m\ A use of Reidemeister-Franz torsion in the combinatorial context also appears in the 
work of Catanzaro, Chernyak, Klein [5] [6]. 

In section TT21 the example of the standard simplex is considered, using and extending 
results of Kalai [2] • 

This paper emphasizes combinatorial identities. We note, however, that as the coefficients 
in the formulas produced are generally shown to be positive integers, a wealth of interesting 
inequalities, often involving simplified expressions, result as unstated corollaries. 

In a future paper we will see how the Cauchy-Binet Theorem can be extended from 
the context of the determinant, which can be regarded as the trace of a one dimensional 
irreducible representation to the context of the trace of an arbitrary irreducible representation 
[32]. By these means, additional combinatorial-geometric formulae can be derived for graphs 
and cell complexes. 


2. Theorem on mesh matrix for d Cycles 

Let X be a finite CW-complex. Orient the cells. This data provides a basis for the 
integral d-chains, Cd(X;Z), a free abelian group. This basis is called the cellular basis. 
Let the boundary map be dj : Cd(X\Z) —> C d -i(X] Z), and set the d-cycles Z d (X\Z) = 
kernel of dd, and ( d — 1) boundaries B d - i(X; Z ) = image of dd, and correspondingly for 
real coefficients R. 

Define a natural integral valued pairing on Cd{X\ Z) by 

C d (X;Z)xC d (X;Z)^Z 

( Zj cij \af , D/j. ^[(T/j,] ) i ^ Zij aj • bj 

where Sd = {[oy]} i s the cellular basis in dimension d. That is, the elements of the cellular 
basis are of length one and orthogonal. 

As above, let rid — |*Sd| be the number of d-cells and Zd, bd be the ranks of the free abelian 
groups, Zd(X] Z), Bd(X ; Z ), of integral d cycles, d boundaries respectively. 
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Introduce the mesh pairings on d cycles and d boundaries: 

/J : Z d {X ; Z) x Z d (X- Z ) C C^A; Z) x C d (X; Z) < A > Z 
z/ : Z) x B d (X- Z) C C^A; Z) x C d (X] Z) < '^>' > Z 

Let d > 1 and £x an ordered basis, say | r = 1, • • ■ ,z d } for the integral d cycles 
Z d (X ; Z) and define the associated mesh matrix for d cycles M{z r } a z d x z d integral valued 
matrix by 

M{z r }i,j =< Zt,Zj > . 

Similarly, fix an ordered basis, say {6 S | s — 1, • • • , b d } for the integral d boundaries B d (X; Z) 
and define the associated mesh matrix for d boundaries M{b z }, a b d xb d integral valued matrix 
by 

M{b s }ij =< bi,bj > . 

In the case of X = T, d = 1, a connected graph, the study of the mesh matrix for 1 cycles, 
i.e., circuits, was initiated by Trent |29j . who showed that in this case: 

det(M{z r }) = ff of spanning/trees of the graph T 

Of course, the determinant of M{z r } is the product of the eigenvalues of the real symmetric 
matrix M{z r }. 

Here Trent’s analysis of the determinant is extended to each dimension of any finite CW- 
complex and to the coefficients of the complete characteristic polynomial, det(t Id — M{z r }) 
of the mesh matrix for d cycles M{z r }. Also, parallel results for M{b r }, the mesh matrix 
for d boundaries, will also be obtained. 

Note that under change of basis of the d-cycles Zj i —> Z k Aj k z k with A G GL(z d , Z), the 
mesh matrix M{z r } changes to the mesh matrix to A-M{z r }-A t , with thus det(t Id—M{z r }) 
changing to det(t Id — A ■ M{z r } ■ A*). ( It is only the determinant, det(M{z r }) which 
is left invariant under all transformations M{z r } i—>• A ■ M{z r } ■ ALThis last holds since 
A G GL(z d , Z) implies that det(A) = ±1.) So it is remarkable that one can arrive at uniform 
geometric formulas for coefficients of the characteristic polynomial recording the dependence 
on the choice of integral basis used. This same observation also applies to det(t Id — M{b r }), 
the characteristic matrix for the mesh matrix for d boundaries. 

Both mesh matrices have a useful decomposition: Let A{z r },A{b s }, respectively, denote 
the integral valued n d x z d , n d x b d matrices representing the inclusions 

i z : Z d (X- Z) C C d (X; Z) and i b : B d (X ; Z) C C d (X- Z) 

with respect to the chosen integral cycle basis {z r } of Z d (X] Z) and the chosen integral 
boundary basis {6 S } of B d (X] Z) and the cellular basis of C d {X\Z). 

Then the respective square z r x z r mesh matrix M{z r } and the b d x b d matrix M{b r } are 
given by: 

M{z r } = A{z r Y ■ A{z r } and M{b s } = A{b s Y ■ A{b s } 

More intrinsically, endow C d (X ; R ), Z^(A"; R ) , B d (X] Z) with the inner products for which 
the specified basis elements, cellular, d cycle, d boundary, respectively, have length one and 
are orthogonal. Denote by . > #z , < .,. > #i> the respective inner products. 
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Let i* : C d (X ; R) —> Z d (X ; R ) the adjoint of i z , and i b : C d (X; R) —> B d (X ; i?) the adjoint 
of i b be defined by : 

< i z (a),b >=< a,i*(b ) > #2 for a G Z d (X;R), b G C d (X]R) 

< i b (a),b >=< a,i b (b) ># 6 for a G B d (X] R), b^C d (X]R) 

Then the entries of M{z r } are just 

AL {z r y U ’ V < f 2 (^u) J f Z ^ Zui'i'z 

showing that the mesh matrix for d cycles M{z r } is just the matrix representing the linear 
mapping: i* ■ i z : Z d (X]R) —> Z d (X]R) and is here given in terms of A{z r } as described 

above. Similarly, the mesh matrix for d boundaries M{b r } is just the matrix representing 

the linear mapping: i b ■ i b \ B d (X]R) —>■ B d (X;R) and is thus given in terms of A{b s } as 
described above. 

Also, one may use the cellular basis, or equivalently the inner product < ., >, to define 
the adjoint d d : C d -i(X; Z) —> C d {X : Z) of the boundary mapping on integral d chains: 

d d :C d {X-Z)^C d _ l {X-Z) 

If K represent the boundary mapping in the cellular basis, the combinatorial Laplacian 
acting on (d-l)-chains is 

= d d • <9* : C d ^{X- Z) C^X- Z). 

It is represented by the matrix K ■ K*. Hence, this setting has a parallel decomposition to 
that for the mesh matrices for cycles and boundaries. This differs slightly from the standard 
usage of this term in differential geometry, as explained in the next paragraph. 

Note: Logically following differential geometry, one could for d — 1 > 0 follow the con¬ 
vention that the Laplacian on d — 1 chains is = d d ■ d* d + where d d _i : 

C d -i(X] R) —> C d - 2 (X] R). However, in view of the combinatorial Hodge decomposi¬ 
tion, Cd_i(A"; R) = R d _i © {image d d } © {image <9qLi}> an orthogonal direct sum de¬ 
composition, the mapping d d maps the eigenvectors of (d d d d \{image d d }) to eigenvectors of 
(' d d d d \{image d d }) with the same non-zero eigenvalues. Since the harmonics R d are identified 
with {c G C d (X;R) | d d (c ) = 0,<9^(c) = 0} and under the combinatorial deRham isomor¬ 
phism map isomorphically to the d th homology, H d (X; R), to record the non-zero eigenvalues 
of A^_! for all d efficiently it suffices to record the non-zero eigenvalues of what is called here 
the combinatorial Laplacian A^_x = d d -d d acting on (d-l)-chains. This observation was used 
in the work of D. Ray and I. M. Singer |26| on their analytic version of Reidemeister-Franz 
torsion. 

In the context of a simple graph, say T, d — 1, C is just the familiar incidence matrix 
from edges to vertices. Moreover, 

Kirchhoff matrix = Graph Laplacian Ao on 0 — chains — C ■ C 1 — d\d{ = Deg — Adj 

where Adj is the adjacency matrix with Adj v>w = +1 if {n,ta} is an edge and = 0 otherwise, 
while Deg is the diagonal matrix with Deg V:V equal to the number of edges containing the 
vertex v. 

For any subset V C S d of the d cells, form the d dimensional subcomplex of X, called Xy 
as follows: 

AV = X (d-1) U V 
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where X 1 is the d — 1-skeleton of X , the union of all cells of dimension at most d — 1. 
By definition, Xy has dimension d, with \V\ cells of dimension d, and the already oriented 
cells of V provide an integral, real basis for C d (Xy ; Z), C d {Xy\ R) respectively. 

For V C Sd, let V c denote the complement of V in Sd, V c = Sd — V. 

Now consider the boundary mapping 

d d :Cd{X-Z)^C d -i(X-Z) 

with image the free group B d ~ i(X; Z) of rank b d ~i- Suppose b d - i > 0. By adding successively 
one by one d-cclls of S d to the d — 1 skeleton of X, X d_1 , the rank the image under d d may 
be increased one by one until the maximal value b d ~ i is obtained. That is, one may exhibit 
subsets V C Sd with 

\V\ = bd- 1 and C d (X v ,R ) -4 B d _\ [X V ] R) c B d -i(X; R) an isomorphism. 

In agreement with the graph case, one calls such a choice of V C S d with these two properties 
a spanning forest of dimension d of X. The kernel of the composite for such V is just 
Z d (Xy ; R), so alternatively, V C Sd is a (d-dimensional) spanning forest precisely if 

\V\ = b d - 1 and Z d (X v ; R) = 0. 

The computation of the determinant of the mesh matrix in dimension d, det(M{z r }), will 
be as a sum over spanning forests of dimension d in X. 

Correspondingly, given a spanning forest V, for each oriented d-cell, say a E V c , the 
boundary d d cr lies in B d -i(X] R) = B d -i(Xy-,R ) so is of the form d d c for a unique c E 
C d (Xy-,Q), by V a spanning forest of X. In this way, a unique rational d-cyclc z[cr\ E 
Zd(X VU { a j; Q) is defined such that its restriction to a has coefficient +1. In particular, 
this shows that if {cri, cr 2 , - • • , cr k } are disjoint oriented cells of S d in V c [with necessarily 
k < |P C | = n d — b d -1 = z<i\, the union 

W — V U {cri, cr 2 , > • • , a fc }, 0 < k < z d 

has the two properties 

|W| = b d - 1 + k and Z d (X w ; R) = k. 

One calls such a set IF C S c j a k augmented spanning forest since it arises exactly by 
adding k d-cells to some (d-dimensional) spanning forest V. In particular, a spanning forest 
is just a 0 augmented spanning forest in this notation. 

Recall of the short exact sequence, 

0 Z d (X ; R) C C d (X; R) % B d _ i(X; R) -E 0. 

Again, in view of this short exact sequence, an equivalent condition for IF to be a k aug¬ 
mented spanning forest is that 0 < k < z d and 

|IF C | =z d -k and j w : Z d (X ; R) C C d {X\ R) re “ C d (X wc ; R) is onto. 

For example, if IF is a k augmented spanning forest of dimension d, then |1F C | = n d — 
(b d - 1 + k) = z d — k. Since the mapping from Z d (X] R) to C d (X W c ; R) has kernel Z d (X] R) D 
C d (X W ] R) = Z d (X w ; R) of dimension k, the image of Z d (X] R) in C d (X W c] R) has dimension 
z d — k. Since |1F C | = z d — k, this implies that Z d (X ; R) —> C d (X W c] R) is onto, proving that 
the first definition implies the second. The converse follows in a similar manner. 
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The computation of the coefficients of the characteristic polynomial of the mesh matrix 
for d cycles of M{z r }, will be in terms of a sum over (d-dimensional) k augmented spanning 
forests with suitable positive integral weights. 

For a finitely generated abelian group G let [G] denote its torsion subgroup and | [G] | the 
order of this subgroup. 

For a finite cell complex X, let t k (X) denote the order of the torsion subgroup of the k th 
homology H k (X; Z ), i.e. , t k (X) = \[H k (X; Z)]\. 


By dehnition, for W a k augmented spanning forest, the mapping 

jw ■ Z d (X ; Z) C C d (X ; Z) -)■ C d (X, X w ; Z) = C d (X W c ; Z) 

is rationally onto. Now using the chosen cycle basis for {z r } for Z d (X ; Z) and the cellular 
basis for C d (Xw y Z), the adjoint is a well defined injection 

f w :C d (X W c-Z)^Z d (X-Z). 

More formally, take the cellular metric < .,. > restricted to C d (Xw c ', Z) and the cycle metric 
< .,. > #z on Z d (X] Z) and form the adjoint. Hence, if jw is represented in these basis by 
the matrix, say Uw , an integer valued (z d — k) x z d matrix, then the adjoint is represented 
by Uf v , the transpose. 

The determinant det(jw • jw) w iH be called the covolume squared of this inclusion jjy of 
the lattice. It may be seen to be equal to the square of the unit volume cell of the sublattice 
inside the real vector space Z d (X; R ) under the metric < .,. > %z . In this case, as Uw is an 
integral valued matrix, this covolume squared is a non-negative integer. 

Similarly, one has the injection of d cycles : 

mw : Z d {X w ; Z) C Z d (X ; Z) 

By W a k augmented spanning forest the free abelian group Z d (Xw', Z) has rank k. Pick 
an integral basis, say of Z d (Xw', Z). For this d cycle basis and the already chosen 

standard basis on {z r } for Z d (X ; Z), let B[W] be the z d xk matrix representing this inclusion 
with respect to the integral basis already chosen. As seen the determinant 

det^BlW]*B[W]) is a positive integer 

and is independent of the choice of integral basis for Z d (Xw', Z). It explicitly records the 
dependence on the fixed basis of Z d (X;Z). 

This determinant is employed in the first main theorem. 

Theorem 2.1. The constant coefficient of (— l) Zd det(t Id — M{z r }) is 

det(M{z r }) = product of eigenvalues of M{z r } = a Zd (M{z r }) 

_ y ( td-i( x V ) \2 

^V, a spanning forest of dimension d \ / 

Also, for V, any spanning forest, i.e, a 0 augmented spanning forest, 

covolume squared of( jf : C d (X V c ; Z) Z d (X ; Z) ) = det(j v ■ jf) = • 

t d - ) 
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If 0 < k < Zd, then 


coefficient of (—l) Zd k t k in det(t Id — M{z r }) 

= (zd — k) th elementary symmetric function in eigenvalues of M{z r } = cr Zd _k(M{z r }) 

= ^W, a k augmented spanning forest C det (B [Wf B [W]) . 


Here the integral matrix B[W] records, as above, the injection mw '■ Zd{Xw] Z) C Zd{X ; Z ). 
Moreover, for W, a k augmented spanning forest, with 0 < k < Zd 


covolume squared of( ;jj v : C d (X W c- Z ) Z d (X ; Z) ) = det(j w ■ j^) 

= ■ dd(B[W\‘B[W\) . 

Additionally, the quotient td 1 ~ 1 ^x) a positive integer for W any k-augmented spanning 
forest; in particular, for k = 0. 


The first part extends Trent’s theorem to cell complexes; the second extends it to evaluate 
the complete characteristic polynomial of the mesh matrix M{z r }. Note that two separate 
expressions are given for the integer weights, one involving X w and the other X W c. The 
explicit dependence on the choice of the integer basis { z r } made for Z d (X]Z) is recorded 
directly through det(B[WYB[W]). 

In fjTj for a geometrically defined basis of Z d (X',Q), given by a choice of d dimensional 
spanning forest,Vo, the coefficients of the complete characteristic polynomial of the cycle 
mesh matrix associated to that geometrical choice are given combinatorially. Similarly, in 
for a (d+1) dimensional spanning forest Vi, there are geometrically defined bases for 
B d (X-,Q). Again in that section, the complete characteristic polynomial of the associated 
geometrically defined cycle and boundary mesh matrices are evaluated combinatorially. 


3. Proof of the main theorem for mesh matrix for d cycles: Theorem 12.11 

In analyzing M{z r } = A{z r }* ■ A{z r } and its characteristic polynomial the following 
elementary standard lemmas are helpful. 

Let A be a m x n matrix. For each pair /, J of subsets / C {1, 2, • • • , m}, J C {1, 2, • • • , n} 
let Ap j be the |/| x \J\ minor of A with row indices from the subset / and column indices 
from the subset J. 

Lemma 3.1. Let A be square n x n matrix, then Ofc(A) the k th , elementary symmetric 
function of the eigenvalues, say Aj, i — 1, • • ■ , n, counted with multiplicities, is given by the 
sum of determinants 

,n} with \I\=k det(Ajy') 

for 0 < k < n. 

Lemma 3.2 (via Cauchy-Binet theorem.). Let A, B be m x n,n x m matrices respectively, 
then the non-vanishing eigenvalues of A ■ B, a m x m matrix, and B ■ A, a n x n matrix 
counted with multiplicities are identical and the k th elementary symmetric function of the 
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eigenvalues of AB or of BA equals: 

@k(yAB^) ... with\i\=k det((AB)ij) 

an d Jc{i, ••• ,n} with \i\=k,\j\=k det(Aj ydetl^B j j') Cauchy Binet 

,m} and Jc{l.-,n} with \I\=k,\J\=k det^Bj^detyAjj) 

= £ Jc with\j\=k det((AB)j ) Cauchy - Binet 
= a k (AB ) 

To prove the first lemma is elementary. The desired coefficient may be evaluated by 
differentiating det(t Id — A) a total of (n — k ) times and setting t — 0. Due to the special 
form of the square matrix (t Id — A), with t's only down the diagonal, this gives a sum of 
determinants of minors along the diagonal blocks of size (k x k) — |/| x |/| indexed by the 
sets (/, I). 

The second lemma follow from the first by utilizing, as in the statement of the theorem, 
the Cauchy-Binet theorem. 

These lemmas can be applied to the mesh matrix on d cycles M{z r } = A{z r Y ■ A{z r } for 
d cycles where A{z r }, a n d x z d integer valued matrix, represents the inclusion: 

i x :Z d (X-Z)cC d (X-,Z). 

The rows of A{z r } are indexed by d-cclls in S d and the columns by indices in {1, • • • ,z d } 
specifying elements of the cycle basis {z r } of Z d (X] Z). 

By the second lemma the (z d — k) th elementary symmetric function in the eigenvalues 
of M{z r } = A{z r Y ■ A{z r } equals the (z d — k) th elementary symmetric function in the 
eigenvalues of A{z r } • A{z r y. By the first lemma this equals the sum of determinants, 

with \I\=Zd~k det( ( A{z r } • A{z r y) {IJ) ), 

where the sum, as indicated, is over the d-cells S d of X which label the rows of A{z r y 
Here there is the equalities of matrices, 

(H{T r } • A{z r } )(/,/) = (^4{- z r})(/,{i,---,z d }) ■ (A{z r } )({i,..., 2 d },j) 

= (^4{-Zr})(/,{l,...,2 d }) ■ ( {A{ z r})(I,{l,--- ,z d }) ) ) 

where (A{z r })^i : {i t ... jZd y' ) represents 

n : Z d (X ; Z) C C d (X; Z) re % ict C d {Xf, Z) 

and its transpose the adjoint for the cycle basis for Z d (X ; Z) and the cellular basis for 
C d (Xf] Z). In toto, this proves: 

The ( z d — k) th elementary symmetric function in the eigenvalues of M{z r } 
with |/|=z d — k det(r j r j j . 


Now clearly, the determinant det(r d ■ rf) equals zero unless 

n : Z d (X; Z) C C d (X; Z) re “ C d (X r , Z) is onto. 

Here |/| = z d — k, so this is precisely the condition that / = W c with W a k augmented 
spanning forest. In this case, 

n = j w : Z d (X ; Z) C C d (X- Z) re “ C d (X we ; Z) . 
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That is, this proves: 

h e (z d — k) th elementary symmetric function in the eigenvalues of M{z r } 

^WcSd a k augmented spanning forest det(j w ■j-fo) . 


Here the adjoint of jw called here j{ v is formed by declaring that the fixed chosen ba¬ 
sis of Z d (X]Z) is orthonormal and that the d-dimensional cellular basis of C d (Xw c ', Z) is 
orthonormal. 


To complete the proof of the theorem for the k = 0 case, for a spanning forest, say V, one 
needs to show that the covolume squared is expressed as 

Now consider 

Jv :Z d (X-Z)^C d (X vc ;Z). 

It is a mapping of free abelian groups and is rationally an isomorphism by the definition of 
spanning forest, so det(jy • jy) is just the square of the order of the torsion abelian group 
C d (X V c;Z)/j v (Z d (X-Z)). ' 

Now quite generally, let IT be a k-augmented spanning tree and consider the restriction 
mapping 

jw ■ Z d (X\ Z) —> C d {X W c- Z) 

which is rationally onto by IT an k-augmented spanning forest, it is claimed that the torsion 
group C d (X W c ; Z)/j v (Z d (X ; Z)) has order 

| C d (X w .| Z)/j w (Z d (X; Z)) I = Vwr , 

M-lV ) 

in particular is an integer as asserted in theorem 12.11 

Hence, in the case k — 0, this implies the desired formula above. 


To see this last write jw as the composite 

jw : Z d (X ; Z) C C d (X- Z) C d (X W c- Z) = C d {X , AV; Z) = Z d (X, X w ; Z) 

since X and X w have the same (d-1) skeleton. Also note that B d (X; Z) maps onto B d (X, Xw ; Z), 
so the desired quotient is isomorphic to the quotient 

H d (X,X w - Z)/ {image H d (X- Z)} . 

But recall from Lyons [18] the utility of introducing the long exact sequence for the pair 
(A", Xw) which reads 

H d {X ; Z) A H d {X, AV; Z) 4 H^Xwl Z) 4 H d _ i(A; Z) ^ H d _ i(A, A w ; Z) = 0 

with a is rationally onto by IT a k augmented spanning forest. 

Hence, letting K denote the image of fj in H d _i(X w ; Z), by exactness, the torsion group 
C d (X,X w ; Z)/j v (Z d (X; Z)) = H d (X,X W ; Z)/ {image H d (X- Z)} is isomorphic to 
image of j3 = K. 

Now consider the associated mappings of torsion subgroups: 

image (5 = K = [K] C \H d _ i(AV; Z)} 4 [H^X- Z)] -> 0 
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By definition 7 ' • j3 = 0. Now let x be a torsion class in H d _i(X ; Z ), say Nx = 0, with 
N a positive integer. By exactness of the pair, there is y G H d _i(X W ] Z) with 7 (y) = x. 
Bnt by 7 (Ny) = 0, Ny is in the image of /3 which has only torsion image. In particular, 
Ny is torsion and so also is y. Hence, 7 ' is onto. Also, this shows the sequence is exact at 
[Hd- i(AV; Z)\. In toto, this proves the desired equality 

\K\ = {[Hd-^Xw; Z)]\/\[H d _i(X] Z)]\ = = a positive integer . 

td- if A J 

The remaining parts of theorem 12.11 will be addressed in the next section 1J4J 

4. Alternative version of the squared covolume det(jw • jw)i cases k > 0 

To complete the proof of theorem 12.11 one must prove 

squared covolume det,(jw ' 3w) = ( d —r^r ) 2 det^B^V] 1 B\W]) 

td 1 (A J 

for W a k-augmented spanning forest, 0 < k < z d . 

To prove this result the following algebraic lemma will play a crucial part. It is proved at 
the end of this section. 

Lemma 4.1. Let 1 < u < n and A be an invertible n by n matrix. Let R be the uxn matrix 
obtained from A by deleting the last n — u rows of A. Let S be the matrix obtained from 
(A -1 )*, the transpose inverse of A by deleting the first u rows, then 

det{R ■ R f ) = det(A) 2 ■ det{S ■ S l ) . 

To be more explicit, as before let { z r } be the fixed chosen basis for Z d (X ; Z), and 

i x : Z d (X ; Z) C C d (X; Z) 

be the inclusion. Use the cellular basis, {1 • a 3 \ j — 1, ■ • • , n d } of oriented d cells to define 
the n d x z d matrix A by 

^-‘crtGSd Ai >r CTt . 

{Here, as is conventional, ix(Z r L r z r ) = E/ c (E r Afc ir A r )o 7 c .} 

Now let W be a k-augmented spanning forest, say W — VUl, with V a choice of spanning 
tree and / C S d — V of size k. 

By V a spanning tree, there is an isomorphism 

j v : Z d (X-Q) ^ C d (X V c-,Q) 

of rational vector spaces. 

Indeed, for each cell cr,- G V c there is a unique rational d-cycle, say z(crj) for X which 
maps under jy to precisely the d-chain 1 • a 3 . That is, 

z(<jj) G Z d (Xy Uc7j ; Q ) restricts to 1 • crj 

for each aj G V c . This geometrically defined cycle z(aj),j — 1, • ■ • ,z d , with aj ranging over 
V c , is uniquely determined by these requirements. 

In this manner, having chosen the spanning forest V, one obtains a canonical rational 
geometric basis { z{<jj ) | aj E V c = S d — V} for Z d (X-, Q ) which restricts under jy to precisely 
the cellular orthonormal basis of C d (Xyc-,Q ). Very nicely the associated natural geometric 
basis for Z d (Xw ; Q) with the same spanning forest V C W is precisely {z(aj) \ <jj G W — V} 
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and similarly Z d (X W c] Q ) has rational basis {z(aj) \ aj G W c }. In this manner, this geometric 
basis respects the direct sum decomposition of Z d (X;Q ): [W = V U I,S d — W = V c U I 
\I\ = k .] 

Q) = Z d (X w ; Q ) © Z d (X VuW cy, Q ) = Z d (X VuI ; Q) © Z rf (AVuvF c )i Q) 

In particular the methods of lemma 14.11 apply for this choice of geometrical basis, those for 
Oj G / and those with oq G IU C . 

As has been demonstrated above in terms of the restriction mapping jx '■ Z d {X ; Z) —$■ 
C d (Xyc] Z) which is rationally an isomorphism with respect to an integral basis of Z d (X; Z) 

i r y \ 

to the cellular one has determinant ± Alternatively expressed, the change of basis 

from the an integral basis of Z d (X ; Z) to the geometrical basis given by the spanning forest 
V has determinant equal to observed above, the k = 0 case of theorem 12.11 

det (A) = ±(^^). 

In a like manner, since V as a spanning forest for Xw, the change of basis from the an 
integral basis of Z d (X\y] Z) to the geometrical basis given by the spanning forest V in Xw 
for Z d (X] Q ) has determinant equal to =k • 


Under this identification of Z d (X] Q ) with C d (X V c ; Q), the mapping j v becomes identified 
with the identity mapping 

Id-.Z d {X-Q)^Z d {X-Q) 

with the the matrix A becoming just the translation between the two bases {z r } and the 
geometric basis z{(Jj) as in 

Z r Ui <j<zj -Aj,r z((J . 

[As is seen by applying the isomorphism jy to both sides.] 

Now arrange that the first z d —k cells of V c are those of W c = S d —W = S d —(V U I) = V c —I 
a total of z d — k in all. Let R be the (z d — k) x z d matrix obtained from A by deleting the 
last k rows. This matrix represents exactly the mapping 

jw '■ Z d (X]Q) —> C d (X W c-Q ) = Z d (Xy uW c-, Q) 

sending z r n- Z,i<j< Zd _ k Aj ir z(<jj ) [ as seen by applying j w \ and det(jw • jw) i s represented 
by the matrix RR T . 

On the other hand, inverting the matrix A, gives 

^(©) Ur (A ^)r,t ^r • 

That is, the identify mapping 

Id:Z d {X-Q)^Z d (X-Q) 

which for the geometric basis on the domain, z((Tj) with G V c and the fixed chosen basis 
{z r } on the range is represented by the matrix A -1 . 

Hence, the inclusion of d-cycle spaces over the rationals 

i\v ■ Z d (X w ; Q ) C Z d (X; Q ) 

is represented for the domain by the geometrical basis z(crj), 1 < j < k a 3 G /, and range 
basis { z r } by the matrix, say S t , and is obtained from A -1 by deleting the first k columns. 
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Equivalently, i* w is represented by S which is obtained form (A -1 )* by deleting the first k 
rows. 

Since H(W] by definition is the matrix for the inclusion Zd(X\y] Z) C Z r j(X ; Z) with 
respect to the an integral basis for Zd(Xw; Z) and the chosen basis { z r } for Zd(X] Z), this 
observation shows that 

det{B[W] 1 ■ B[W ]) = ( ? d -f* v) . ) 2 ■ det(i* w ■ i w ) 

td—a v^- W ) 

and iw • i\y represented by 

( )r,t\l<t<Zd—k ) ■ (A )r',t'\l<t'<z d —k 

= ( ((dl- 1 )^ U< t < Zd . k ) • (((A~ 1 ) t ) t t ,y | l<f<z d -kY 

= s • s*. 

By lemma 14.11 

det(R ■ R f ) = det(A) 2 ■ det(S ■ S l ) . 

Hence, one sees that, in the notation of theorem 12. II 
det{ iw ■ jw ) = det{R ■ R f ) = ( det(A )) 2 ■ det(S ■ S*) 

= (■VTirf) 2 ' (svSr) -2 ' det(B\W]‘B\W]) = ■ dei(S[lT]‘fl[lT]) 

As noted above, H[hE] is an integral valued matrix and is an integer when W is 

a k augmented spanning forest. So this completes the proof of theorem 12.11 

Proof of lemma 14.11 

One first notes that by the Cauchy-Binet theorem the determinant of R- R t equals the sum 
of squares of the r x r minors of A which occur in the first r rows. Similarly, the determinant 
of S ■ S f equals the sum of squares of the s x s minors of (A -1 )* which occur in the last s 
rows. Hence, the above lemma will follow from the following lemma: 

Lemma 4.2. Let A be a square (r + s) x (r + s) non-singular matrix, i.e, det(A) ^ 0. Let 
E be a r x r minor of A and F be the complementary s x s minor of then 

det(E) = ±det(A ) • det(F) . 

Now replacing A by pre and post multiplication by permutation matrices, with care about 
signs, this last lemma is reduced to: 

Lemma 4.3. Let A be a square (r + s) x (r + s ) non-singular matrix, i.e, det(A) ^ 0. Let 
E be the top left r x r block of A and F be the bottom s x s block of (A*) -1 . Then 

det(E ) = det(A ) • det(F) . 

By continuity, it suffices to consider the case with det(E) ^ 0, det(F) ^ 0. For these cases 
this lemma follows from the standard theorem. 

Theorem 4.4. Let A be a square (r + s) x (r + s) non-singular matrix, of the block partitioned 
form: 
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where a, b, c, c are therxr, rxs, sxr, sxs matrices respectively, defined by a = ( Aj|i<i<r, 1 <j< r 

( r+l<j<r+s );C ( |r+l<i<r+s, l<j<r |r+l<i<r+s, r+l<j<r+s ) respec¬ 

tively. 

If det(a) 0, then det(A) = det(a) ■det(d — c a~ 1 b ). In particular, the inverse (d —c a -1 b ) _1 
exists in the context det(A) 0 ,det(a) 0. 

Ifdet(d) 0, then det(A) = det(d)-det(a—b d _1 c ) In particular, the inverse ( a—b d -1 c) _1 
exists in the context det(A) 0, det(d) 0. 

If det(A) 0 ,det(a) 0 ,det(d) 0, then the inverse of A, A~ l , has the block decompo¬ 
sition : 


A -1 


(a — b ct _ 1 c ) _1 
— (d — c a^ 1 b ) _1 ca~ l 
( (a — b d _ 1 c ) _1 

l —d _1 c ( a — b d~ l c)~ l 


— (a — b d^ 1 c) -1 bd~ l \ 

(d — c a -1 6) _1 J 
—a~ l b (d — c a^ 1 6) _1 \ 
(d — c a~ l 6) _1 J 


),b 


The proof of this is elementary. By det(a ) ^ 0, one may reduce A to Schure complementary 
form by multiplying by a matrix of determinant 1; as in: 


fa 0 \ _ f a b \ f Id —a 1 b \ 

y c d — c a -1 b J y c d y y 0 Id J 

Hence, det(A) = det(a ) • det(d — c a -1 c b ) as claimed. Similarly for the case det(d ) ^ 0. 
The formulas for A follow by direct calculation of A~~ l ■ A. 


5. The mesh matrix for d boundaries; evaluation of its characteristic polynomial 

There is a parallel development for the mesh matrix on d boundaries. 

One fixes a choice of basis for the integral d-dimensional boundaries B d (X;Z) say {6 r } 
and declares that they form an orthonormal basis for the boundaries. 

Let V C S d be called a spanning coforest if the restriction mapping 

k v : B d (X ; Q) C C d (X; Q ) re “ C d {X v \ Q ) 

is an isomorphism. In particular, \V\ = b d = rank B d (X ; Q). This is called by Lyons [18] a 
cobase. 

By incrementally making appropriate choices of d-cells in S d , the image of the restriction 
can be made larger and larger until it achieves the maximal value b d . Hence, such spanning 
coforests always exist. 

Now let W be the result of deleting k cells [ with 0 < k < b d ], from a spanning coforest, say 
V. The result is called a k reduced spanning coforest. It is characterized by |W| = b d — k 
and 

k w : B d (X ; Q) C C d {X- Q ) re “ C d {X w - Q ) 

is surjective. Naturally, a spanning coforest is just a 0 reduced spanning coforest in this 
terminology. 

This k\y is obtained by tensoring with the rationals the composite mapping 
k w : B d (X ; Z) C C d (X; Z) ->■ C d (X, X Wc -Z) = C d (X w ; Z), 
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which is rationally a surjection by W a k reduced spanning coforest. Its adjoint with respect 
to the orthonormal basis {b r } of B r (X\ Q ) and the cellular basis [declared to be orthonormal] 
of Cd(X\yc] Z ) is the inclusion : 

k* w :C d {X w -Z)^B d {X-Z). 

Again, one calls the determinant det(k\y • kyy) the covolume squared of this inclusion kf v 
of lattices. 


Analogously to the mesh cycle cases, let {6 S } be a fixed choice of integral basis for B d (X ; Z) 
and introduce the boundary mesh matrix 

M{b a } = {<b i ,b j >} 

where one includes B d (X] Z) C C d (X;Z) and used the cellular basis to define the inner 
product on d chains. Equivalently, this boundary mesh matrix is kx ■ k* x . with the adjoint 
computed with respect to the respective basis {6 S } on B d (X;Z) and the cellular basis for 
the d chains made orthonormal in the respective inner products. 

As seen before, the application of lemma 13.21 evaluates the b d — k elementary symmetric 
function of the eigenvalues of the boundary mesh matrix M({b s }) as a sum : 

a b d -k(M({b s })) = W ith \w\=b d -k det(k w ■ k^) 

Again the determinant is non-zero only for the case that kyy is onto, that is the case that W 
is a k reduced spanning coforest. Hence, 

k (Af ({hs }) ) £lV with W a k reduced spanning coforest dd( k\y • kyy^j . 


For a k-reduced spanning coforest W chose a spanning coforest V with W C V. That is, 
W is obtained from V by deleting k oriented cells, say I,soV — W U I and \I\ = k. Having 
made such a choice of V, consider the inclusion mapping : 

X\v '■ B d (X W c] Z) C B d (X ; Z) 

For a choice of integral basis for B d (Xw^', Z) = B d {Xv^r, Z)) and the chosen basis {b s } 
for B d (X]Z) let H'JIF] denote the associated integral matrix. Hence, the b d x b d matrix 
B'lWYB'lW] represents x^y-x^ and det^B'lWYB'[W]) is independent of the choice of integral 
basis of B d (X W cZ) = B d {Xyc { ji\ Z)). 

By W a k reduced spanning coforest equal to lb — V — I,\I\ = k for a spanning coforest 
V, the restriction mapping 

Uw '■ B d (X W c ; Z ) = B d (X V c Ul ] Z) (->• C d {Xf, Z) 

will be a rational isomorphism. 

The parallel result to theorem 12.II for the mesh matrix for boundaries M{b s } is: 

Theorem 5.1. For V a d dimensional spanning coforest of X, the group H d (X; X V c] Z) is 
a torsion subgroup. Let its order be denoted by v(V,X). 
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The constant coefficient of (— l) bd det(t Id — M{b r }) is 

det(M{b s }) = product of eigenvalues of M{b s } = a bd (M{b s }) 

a spanning coforest of dimension d (u(P,-Xj) 

Also, for V, any spanning coforest, i.e, a 0 reduced spanning co forest, 

covolume squared of( kf : C d (Xy] Z ) B d (X: : Z ) ) = detfky ■ kf) = (v(V,X)) 2 . 


Quite generally, let W be a d-dimensional k reduced spanning coforest in X. For such W, 
the quotient C d (X W c] Z) / Xw(Z d (X\ Z)) is a torsion group isomorphic to H d (X-,X W c] Z). 
Let the order of this finite torsion group be denoted by v(V,X). 

Moreover for W contained in a spanning coforest V with W = V — I for |/| = k, then the 
quotient C d (Xj; Z)/ y w (B d (X W c)-, Z)) is a torsion group and is isomorphic to H d (X W c , X V c ; Z). 
Let the order of this finite torsion group be denoted by u(W, V). 

It is claimed that the quotient 


fiW) = 


u(w .; v) 


v(V,X) 

is independent of the choice of V and moreover f(W) v(W, X) is an integer. 
If 0 < k < b d , then 


coefficient of (— l) bd k t k in det(t Id — M{b r }) 

= (b d — k) th elementary symmetric function in eigenvalues of M{b r } = a bd _ k (M{b r }) 


= E 


W, a k reduced spanning coforest 


'u(w,v )n 

. v{V,X) 


det((B'[W]yB'[W}) 


Moreover, for W, a k augmented spanning forest, with 0 > k < z d : 


covolume squared of( kf v : C d (X w \ Z) c —t B d (X ; Z) ) 

= (^) 2 det((B’{W])‘B'{W]) 

In particular, the coefficients of the above sum over k reduced spanning coforest yielding 
a bd -k(M{b r }) are all integers. 


Note that if H d (X;Q ) = 0, then Z d (X;Q) = B d (X;Q ) and each spanning forest W has 
complement exactly a spanning coforest W c and the complement of a k augmented spanning 
fores is a k-reduced spanning coforest. In particular, in this situation theorems 12.11 and 12.11 
become identical. 


6. Proof of the main theorem for mesh matrix for d boundaries: Theorem 15.11 


Now consider the case V a spanning coforest, i.e., k — 0, the mappings ky, ky are rationally 
isomorphisms. Hence, the determinant detfky ■ ky) is the square of the order of the finite 
abelian group: 

C d (X v ; Z)/{ image by k v of B d (X ; Z) } 

It is claimed that this finite torsion group has order v(V,X) completing the proof for the 
case k = 0. 


Quite generally, let W be a k reduced spanning forest, then k w : B d (X ; Z) —> C d (X v \ Z) 
is rationally onto, so the quotient C d (Xw', Z)/{ image by kw of B d (X]Z ) } is a torsion 
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group. It is claimed that the order of this finite quotient group 

| C d (X w ; Z)/{ image by k w of B d (X ; Z) } | 

is v(W,X). 

The mapping kw may be reinterpreted [ since X w , X W c each have the same d — 1 skeleton 
as X : 

k w : B d {X ; Z) ° 4 ° B d (X, X W c ; Z) C Z d (X, X W c ; Z) = C d (X w - Z ) 

So this finite abelian quotient is identified with H d (X, Xwc; Z). Its order is v(W,X). 

Finally, in view of the long exact sequence of the triple (X, Xwc, Xyc), there is an exact 
sequence 

H d (Xw c , Xyc\ Z ) — > H d (X, Xyc ; Z) —)■ H d (X, Xw c ', Z) — > H d ^i(Xw c , Xy c ; Z) = 0 

since Xwc, Xyc have the same (d-1) skeleton. By the above all these are finite torsion groups. 
Hence, by exactness, u(W, V)/(v(V,X)/v(W,X)) = ^ -v(W,X) is an integer as desired. 


Next consider again as above V a spanning coforest with W = V — I for \I\ = k oriented 
cells of V. That is, W is a k reduced spanning coforest. Now the mapping 

Uw '■ B d (X W c'i Z ) = B^Xyctjj] Z) —> C d (Xj ; Z) = C d (X V c Ull Xyc] Z) 


is rationally an isomorphism. Recall that W c = V c U I. It can be factored as 


B d (X W c- Z) ° 4 ° B d (X w c, AVc; Z) C Z d (X w c, X V c ; Z) = C d (X w c, A yc ; Z) . 


Hence, the quotient C d (Xi] Z)/yw{B d (Xw c ', Z)) is a finite torsion group and is isomorphic 
to H d (X W c,Xyc] Z). Its order is u(W, V). 

Let A! be the n d x b d matrix represent the mapping kx '■ B d (X;Z) C d (X;Z) with 
respect to the chosen basis {b s } for B d (X\ Z) and the cellular basis for C d (X; Z). In these 
terms the boundary mesh matrix, M({b s }), is just {A!)* ■ {A!). 


Now to evaluate the covolume det(kw • kw) for W a k reduced spanning coforest with 
0 < k < b d , pick a spanning coforest V of X with W C V. Let W = V — I with |/| = k. 
Then by V a spanning coforest the restriction mapping 

k v : B d {X- Z) C d (Xy-Z) 

is rationally an isomorphism. Also by above, the quotient is of finite order v(V,X). 

In particular, ky is rationally an isomorphism, so for each <jj G V [oriented as in S d ] there 
is a unique d boundary b(aj) G B d (X;Q) which restricts to 1 • a in C d (X W ',Q)- That is, 
b(<jj) lies in B d (X V c Uaj ; Q ) Here <7j ranges over V which is of size b d , so 1 < j < b d . 

Linder the identification above of B d ( X] Q ) with C d (X v ), identifying b(aj) with 1 • a 3 for 
<jj G V, the map ky is identified with the identity mapping 

Id : B d (X] Q) — >■ B d (X; Q) 

with A' recording the change of basis mapping 

b s — ^<Tjev,j=i, -b d . 
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Note that there ia a direct sum decomposition 

B d (X- Q ) = B d (X V c UW , Q)) ® (B d (X- Q) n C d (X V c n7 ; Q) 

which is compatible with the orthonormal basis given by the geometric generators {b(<Jj) \ Oj £ 

V}. The first summand has basis { b(aj ) j = 1 • • - b d — k} = { b(aj ) | a 3 £ W}; the second 
summand has basis {b(aj) \ a j £ I}. 

Arrange that the first b d — k cells of V be those of W = V — /, and let R' be the 
matrix obtained from A' by deleting the last k rows. Then R' represents k\y and k\y ■ k^ is 
represented by R'(R'y. 

As before let S be the last k rows of ((A') -1 ) 4 . Then S 1 represents the mapping 

Bd(Xyy c \ Z) = B d (Xyc n r] Q ) C B c j(X ; Q) 

with the geometric basis used as a basis for the domain and the the fixed chosen basis {fe s } 
for the range. 

Converting from the geometric basis on the domain to any choice of integral basis for 
B d (X\yc ,; Z) entails a change with determinant the same as the order of 

C d {X 1 -Z)ly w [B d {X wa -Z)) 

which is by definition u(W, V). That is, 

detiStf) = u(W,V) 2 det(B[W}'(B[W}'Y) 

Appealing to lemma 14.11 gives the desired result 

det(kw ■ k^y) = det(R' ■ ( R'Y ) 

= det(A') 2 det((S')(S'y) = det(A') 2 u{W, V)~ 2 det{B'[W](B[W]'y) 

= (u(W,V)/v(V,X )) 2 det{B'[W]{B[W]'y) 

as desired. 

Since the determinant det[kw-kyy) as a covolume of an integral lattice in the standard lat¬ 
tice with standard inner product is an integer, the quantity ( u{W , V)/v(V, X)) 2 det(B'[W](B[W]'y) 
is an integer also. 

The equality det(kw • k^) = (u(W, V)/v(V, X)) 2 det(B'[W}(B[W}'y) proves that the quo¬ 
tient f{W) = is independent of the choice of V. 

7. Evaluation of all coefficients of mesh matrices for cycles and boundaries with 

respect to geometrically defined bases. 

Let Vo C S d be a fixed choice of d dimensional spanning forest of A". As explained in §[2l 
there is a natural geometric basis for the rational d —cycles Z d (X;Q) described as follows: 

For each oriented d cell, say <jj in S d — Vo, let z(<jj) be the unique d cycle in Z d (F U o 3 \ Q ) 
which evaluates to +1 on o 3 . As explained there, as o 3 ranges over S d — Vo the d cycles z(a 3 ) 
give a basis for the rational vector space Z d (X;Q) of d cycles on X. Also, for any subset 
U C S d — Vo, the elements {z(o~j) | <Jj £ U} forms a rational basis for Z d (V 0 U U; Q). 

Associated to this geometric basis {z(<jj)} one defines the cycle mesh matrix of A" for 
this geometric basis using Vo as M cyc i e ^ Voj x){z((Tj)} via (a, a') i —>< z(a),z(a') >x with the 
pairing defined using the cell complex A". 
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Theorem 7.1. Let Vo be a d dimensional spanning forest, then the cycle mesh matrix 
McycledL,x){z{Cj )} is a z d x z d matrix with coefficient of {—l) k t Zd ~ k the k th elementary sym¬ 
metric function of its eigenvalues, called here oy(Vo). Moreover, there is a direct combina¬ 
torial formula for this coefficient: 


&k(y o) — £ UcS d -Vo with \U\=k (^V, a spanning forest in Xy 0 


td lAV) N2 
td-l( X V 0 l>u) ' 


Analogously, let V\ C Sd+i be a fixed choice of d + 1 dimensional spanning forest of X. 
Again, there is a natural geometric basis for the d boundaries B d (X]Q). It consists of the 
boundaries {<9d+i(l -u) \ u G Sd +1 — Vi}- Also, for each subcomplex of the form U V\ \JU 
with U C Sd +1 — V , the images {^ + i(l ■ u) \ u G S d +i — V and u G U} form a rational 
basis for B d {X^ U Vi U U\Q), the boundaries. 

Associated to this geometric basis {d d +i(u) \ u G S d +i — V }, one defines the boundary 
mesh matrix of A" for this geometric basis defined by V as 

M boun dary,(y 1 ,x)({dd+i{u) \ u G S d + 1 — Vi}) via (u,v) HG< dd +1 (u),d d +i(v) > x with the pairing 
defined using the cell complex X. 

Theorem 7.2. Let V\ be a d + 1 spanning forest, then the boundary mesh matrix 
M boundary ^ Vu x)({d d +i(u) | u G S d + 1 - Vi}) is a b d x b d matrix with coefficient of (-1 ) k t bd ~ k 
the k th elementary symmetric function of its eigenvalues, called here oy(Vi). Moreover, there 
is a direct combinatorial formula for this coefficient: 

&k(f l) ^ L/C Vi with \U\=k ( ^V, a d dimensional spanning coforest in X^UViLiU I ^d— 1 (-^-Vi t ^-W c i Z') | ) 

where \H d -i(Xv 1 , X\yc; Z)\ is the order of the finite torsion group H d -i{Xv 1 , Xyrc) Z). 
with X Vl = X& U V 

and X W c = A^ d_ b U (S d — W). This weight is alternately expressed as 
the integer [ 1 ^ d ~ 1 ^ Vl ' > ] ; see the next section. 

The proofs of these two theorems proceeds as follows: 

For theorem 17.11 applying lemmas 13.1113.21 one sees that ay (Vo) is expressed as a double 
sum over U C S d — Vo with \U\ — k and over Y C S d of size A; of a square of the determinant 
for the restriction mapping Z d (V 0 U U;Q) to C d (Y ; Q). This last mapping is onto [and hence 
an isomorphism] only if Y c (Vi U U) and V — (Vi — Y) is a spanning forest for X VlUU . 

In this instance, the matrix at issue is that for changing from the geometrical basis for the 
spanning forest Vo in X VoUU to the geometrical basis for the spanning forest V of the same cell 
complex Xy 0 \ju- These two bases can be compared to the integral basis of Z d (Xv 0 uu] Z). by 
invoking theorem 12.11 On that basis, the desired determinant is W ( td -V x v) \ _ 

° - y td-l( X V 0 UU) ’ V td-l(Xv 0 uu) ' 

td-i( x v) 

td-l( x v 0 ) ' 

For theorem 17.21 applying lemmas [371113.21 one sees that oy(Vi) is expressed as a double 
sum over U C V\ with \U\ — k and over V C S d of size k of a square of the determinant for 
the mapping 

C d+ i(A< d) U V h Z) B d {X {d) U V, Z) C C d (X; Z) restr Z? ion c d ( A (d) U V, Z) . 

This last mapping is onto [and hence an isomorphism] only if V C X^ d> u JJ and V is a 
d dimensional spanning coforest for A*^) U U. As the next section proves, in this case the 
mapping is rationally an isomorphism with quotient group a finite order equal to the order 
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of the torsion group H d _i( X Vl ,X W c; Z ), Here X Vl 
This is alternately expressed there as the integer 


= X& U V x and X W c = X^ U(S d -W). 

I *d-l Awkd-2 A(Vi) c kd-iO H 1 

L td- 2 P0 


8 . A Kirchhoff type theorem for cell complexes, extending Lyons’s [ 18 ] work to 

the complete characteristic polynomial 


In the same manner: 

Let A be the integer valued matrix representing the boundary mapping with respect to 
the cellular basis 

d d :C d (X-Z)^C d _ 1 (X-R). 

Then A t represents the adjoint for the cellular inner product < Recall that d d has rank 
b d - i = dim B d _, (X ; R ) 

Let 0 < m < b d - 1 with S d the d cells with their specified orientations. 

By the same algebra as before, one gets 

m th elementary symmetric function of the eigenvalues of d ■ d* 

— ^ics d -1 Jcs d with \i\=\j\= m [ det{A] 2 . 

For the term to be non-zero, the mapping Ayj) must be an isomorphism after tensoring 
with the reals. 

The matrix Ayj) represents the composite : 

C d (Xj ; Z) C C d (X; Z) 4 B d _ i(A; Z) C C^X; Z) re % ict C d - X {.Xy Z) 

For det(Ay j')) 0, necessarily the composite 

C d (Xj- R ) C C d (X- R) 4 B^X- R) 

must be one to one. That is, Z d (Xj ; Z) = 0. 

The conventional notation for this is : J C S d is called a forest of size | J| = m. That 
is, Xj carries no d-cycles. 

Additionally, for det(Ayy)) ^ 0, necessarily the composite 

Ba-^Xj- R) C B d _ x {X- R) C C d - X {X- R) res % ict C' d _ 1 (A J ; R) 

must be an isomorphism. That is, / C S d -\ interpreted via Xj = X d ~ 2 U / inside Xj = 
X d ~ x U J must be a spanning coforest of Xj. In combination this implies that |/| = m also. 

Hence, the above sum is restricted to J = V a (d-dimensional) forest of size m and I = W 
a ((d-l)-dimensional) spanning coforest of Xj = Xy. This ensures that det(Ayy)) ^ 0. 

m th elementary symmetric function of the eigenvalue of d ■ d* 

^ V<ZS d WcS d —i with | V\=m, | W\=m, V a forest of size m and W a spanning coforest of Xy [ det^Ayyy')') ] 

Consider the inclusion X w = X d ~ 2 U W C Xy = A d_1 U V, defining the mapping A^ WV) 

C d (Xy-Z) C C d (x- Z) 4 B d _ x {X- Z) C Cd-^X-Z) res d\ ict C d -i(X w ; Z) 
which in the cases of interest is a rational isomorphism. 
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Now as A. Duval, C. Klivans, and J. Martin note [T2] . one can reinterpret this map more 
geometrically. On integral chains, 

C d {X v , X W c- : Z ) = C d (X v] Z) 4 C d ^{X v ) re “ C d ^( AV; Z) = C d _i( X v , X w <; Z) 

becomes identified with the boundary mapping for the pair (Xy, X'w c )- Since C d -\{{Xy, XZ) 
0 and C d +i((Xv, X\yc' : Z) = 0 and this boundary mapping is a rational isomorphism, it fol¬ 
lows that 

det{A(y W )) 2 = \Hd-i(Xv, Xjyc; Z)\ 2 

in the cases of interest. 


In the maximal case k = b d , by V a forest of size b d , necessarily the mapping C d (X v : Z) C 

C d (X] Z) 4 B d _i(X ; Z) is a rational isomorphism. Hence, the set of W for which W is a 
spanning coforest of Xy is the same as the spanning coforests of A". 

Lyons proved in the case k = b d , that 


|H d _ 1 (AV,AV;4l 2 = [ 


0-i {X w )t d - 2 (. X V c )f^_ 1 (V) 
t d - 2 (X) 


G A 


where for a CW complex K the order of the torsion of Hk(K ; Z) is denoted by tk(K) and 
t'd^iy) is the order of the quotient of Z d _i(X; Z) by Z d _fX) D B d _ i(X; Q) + Z d (X vc ; Z). 
Here by convention Xw = X ( ' d ^ l) U W and Xyc = X^ U V c . 

For m = b d the corresponding sum decouples, yielding the first part of the next theorem, 
which appears in the work of Lyons: 


Theorem 8.1. a. Lyons [T8 j. The product of the non-zero eigenvalues of d ■ d*, the 
combinatorial Laplacian, or Kirchhoff mapping, acting on (d-1) chains C d _i(X\ Z) —» 
C d _i(X] Z), i.e., equals 

— f Y (i (Y H 2 'l l Y n 2 \ 

V ^W , forest of size m=bd y^d— )) J * \ spanning coforest in X \ td—2(X) / /’ 

where each term of each sum is an integer. 

b. More generally, if 0 < k < b d , the k th elementary symmetric function in the eigenval¬ 
ues of d ■ d*, the combinatorial Laplacian, or Kirchhoff map is 


= E 


V a 


forest of size k 


E 


W a spanning coforest of Xy 


H d {Xy,X W c-Z) I 2 ] . 


9 . Homology covolume squared via determinants of mesh matrices 


Naturally, the above results imply a statement about H d (X] R). Recall that using the cellu¬ 
lar basis to identify chains and cochains, the coboundary dual to d d : C d (X ; R ) —> C d -i(X] R) 
gets identified with d d : C d _i(X : Z) — > C d (X]Z). The Hodge theory decomposition be¬ 
comes in this setting an orthogonal direct sum decomposition of real d chains with the cellular 
inner product < ., >: 

C d (X;R ) = Tid © {image d d+ i} © {image d d }, 

with Ti d the space of (d-dimensional) harmonics TL d = {c \ d d (c) = 0 and d d+1 (c) = 0}. 

Now one has the direct sum orthogonal decomposition : 

Z d {X-R) = n d ®B d {X-R) 
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Let B d (X] Z ) = B d (X ; Q) ; Z). Then B d (X ; Z) is a direct summand of C d {X ; Z) and 

the quotient B d (X ; Z)/Bd(X ; Z) is isomorphic to the torsion in Hd(X ; Z). 

Thus the covolume squared for the lattice B d (X ; Z) in C^pf; Z) with cellular inner product 
< ., > is t d (X) 2 times that of B d (X] Z). 

Now chose a direct summand, say G C Zd(X\ Z), with 

Z,(X;Z) = G®R d (X;Z) . 

This is feasible as B d (X ; Z) is a direct summand. Then the covolume of Z d (X] Z) computed 
from this lattice is exactly the same as that of tt(G)®B d (X; Z), where 7T d is is the orthogonal 
projection of Z d (X ; Z) onto H d or equivalently the direct summand G. This shows that 

covolume squared for lattice Z d {X ; Z) under < > 

= ( covolume squared for lattice B d (X ; Z) under <.,.>) 

•( covolume for the lattice n d (Z d (X; Z)) C 40 under <.,>). 

That is, if one defines the lattice in H d by the orthogonal projection n d of the lattice 
Z d (X] Z) onto H d , one has 

covolume squared for lattice TT d (Zd(X-, Z)) C H d wider < .,. > 

_ covolume squared for lattice Zd(C\Z) under 

covolume squared for lattice Bd(C\Z) under / 

This gives in view of the above a combinatorial description of this “homological covolume”. 

10. Combinatorial identities via a computation of Reidemeister-Franz torsion 

As explicated by Milnor [21], the R-torsion of Reidemeister-Franz [27] [13] 
t(C*, {di}, {ejj}, {hj t k}), is an invariant which is a non-zero real number, associated to the 
following data: 

• A finite chain complex of real finite dimensional vectors spaces, say Co, • • • , Cn with 
differentials di : Ci —>■ Cj_i with di-\ ■ di = 0 as recorded by the sequence: 

0 -> C n h C ni ■ ■ ■ 4 C 0 -)• 0 

Let Ci have real dimension r\. Let Z* = {kernel of di : Ci —^ Cj_i} be the i-cyclcs; 
let Bi = {image of d* + i; Cj + i —> C*} be the i-boundaries. Let Hi be the i th homology, 
Hi = Zi/Bj and have rank s t . Let Bi have rank t t . 

• For each vector space Ci a chosen basis, say {e^- | j — 1, • • ■ , rj} 

• For each homology Hi a chosen basis, say | k — 1, • • • , s;} 

This torsion r is defined as follows: Chose elements hi^ € Zj C Cj which project to the 
chosen homology classes Next chose elements, say jA-iq | l = 1, • • • , U-i} which are a 
basis for Bi- 1 . Next chose elements G Ci with di(&ij) = 6,- —ip with l — 1, ■ ■ • 

With these choices, Z* has basis U {fej,;}] and C,; has the basis 

Xi = [{fh.k} U {b it i} U {c iiV }\. Now set D t = \ det([Xi : {e^}]) | where [X t : {e^}] is the 
matrix expressing the basis X t in terms of the basis {eij}. One defines the Reidemeister- 
Franz torsion via : 

N 

(r(C„ {d,}, {e (J }, {ft M })) 2 = n (A)'- 11 * 

i=0 
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A simple check shows that the R-torsion r squared is independent of the choices of h i ^ 1 h- 1 , 1 , c i,v 
[26] ; also, the sign issues are disposed of by taking the square. 

Turning to the case of a cell complex of dimension d, one has the cellular basis for C t (X ; R) 
given by the i-cells of X, once oriented. Moreover, the real homology Hi(X]R) has a pre¬ 
ferred basis coming from a choice of Z-basis, say {h, t k} for the free part of the integral 
homology. Consequently, there is a canonically defined R-torsion for any finite cell complex 
A", t(C'*(A; R),dji, {combinatorial basis for C*(X;i?) and H+(X; R)}). Not surprisingly, 
there is a simple formula (known to experts) for it: 

Theorem 10.1. The combinatorial R-torsion is given by 


d -1 

[ t((7*(A; R) , 8r, {combinatorial basis for C* and H *)} ] 2 = [ W i 2 • 

k =1 

where 9 x ,i is the order of the torsion in HfX] Z ). 

Proof: Since Z t (X ; Z) is a free direct summand of C t (X] Z), the quotient 
Ci(X; Z)/Bi(X\ Z) has torsion subgroup equal to the torsion subgroup of Hf X] Z) = 
ZfX] Z) / BfX; Z)] so in particular its order is just 9 Xji . Alternatively expressed, let 
Bf{X; Z) be the lattice of elements x E Cf X; Z) with Nc E B^X ; Z) for an integer N ^ 0, 
then \Bi(X] Z)/Bi(X\ Z)\ = 9 x ,i and B.f A"; Z) is a direct summand of Ci( A; Z). 

Now chose elements say 6^ E Bf X\ Z) that are a Z basis for the direct summand B^X ; Z). 
Since Zi(X; Z)/B i (X\ Z) = H^X; Z)/ {torsion} is free abelian, one may choose classes 
in Zi(X]Z) which project to a Z-basis for Hi (A; Z)/{ torsion}, i.e., the required 
combinatorial basis. With these choices made {h h k} LI form an integral basis for the 
direct summand Z-lattice ZfX ; Z) inside C*(A; Z). Again, since Cj(A"; Z)/ZfX ; Z) is free 
abelian, one may chose elements Xi ^ E Ci(X;Z) so that they project to an integral basis 
of this quotient. With these choices the lattice CfX]Z) has an integral basis given by 
U {bij} U {xj^}. This gives a reference basis for Cf X] Z) and represents in Cf A"; R ) 
a basis which differs from the cellular basis by an element in GL(dim CfX] R), Z) so the 
transition matrix has determinant ±1. Consequently, one may use this reference basis in 
computing the Dt and the R-torsion r. 

Naturally, one uses the chosen lift h .^ again and the classes {di + i(xi + i^)} and {xi^} in 
computing the R-torsion with the above reference basis. That is, 

Xi = [{hi t k}, {di+i^i+i,;')}, }] is a basis for Cj(A; R) which conforms to the above recipe. 

Xi is a sublattiice of the reference lattice [{hy*,}, {&*,«'}, of index |Rj(A; Z)/Bf A; Z)\ = 

9 Xji . Hence, Di = 9 Xti and the theorem follows. 

Returning to the general context of the R-torsion, a chain complex with the above choices 
{ x i,j}> {hi,k}i there is a real inner product in which these basis elements are orthonormal. 
Hence, for these inner products the boundary mapping d i+ \ has an adjoint 

d* +1 : Ci Ci + 1 . 

As before let {hi^} be a specified basis for the real vector space Ci. 
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As Ray and Singer remark in their foundational paper on analytic torsion [26], there is a 
formula for the R-torsion r in terms of the reduced Laplacian 

A.y ,i ■— (di +1 • d*\B,) : B, — y B,, . 

The first part of this relation is a measure, via covolumes of the chosen homology basis: 
For this purpose, note that one has an orthogonal Hodge decomposition, 

Ci = Harrrii © {image of dj+i} © {image of d*_ ± } , 

in which Harmi = {h \ di(h ) = 0 and d*_ 1 (/r) = 0} are the i-harmonics. Moreover, Z t = 
Harmi © {image of d} is a orthogonal decomposition. Let 7q be the orthogonal projection : 


7 Ti : Ci —)■ Harmi 


It induces the isomorphism HfX ; R ) = Z t /B t Harmi . 

In this manner the specified basis for the integral homology modulo the torsion {/r^*.} 
under the orthogonal projection to Harmi yields the desired elements 7r i(hi^) G Harmi. 
They define a Z-lattice of full rank which are independent of the lifts chosen for the {A^j’s. 
The subspace Harmi inherits an inner product from that on Ci. One defines 

{'Hx,i) 2 = [covolume(lattice { 7 Ti(h i}k )} in Ci )] 2 = det({< TTi(hi, a ),^i(hi,b) >}) 

as a measure of the size of the combinatorial homology basis. 

With these preparations one has: 


Theorem 10.2 (Ray and Singer [26]). The R-torsion is given in terms of the covolumes Hi 
for the i-homologies and the determinant of the reduced Laplacians by 

[ r(C*(X; R ), 3r, {combinatorial basis for C*{X]R) and L/*(A;i?)) ] 2 

= n“,o [*f(A.v,,) ■ ■ 

[Note that in Ray and Singer’s paper they use the full Laplacian A = d l+ \ d* + d t 
with corresponding modification of the above formula. They also restrict to the acyclic cases 
Hi = 0 .] 

Returning to the cell context and the combinatorial bases for chain and homologies, The 
formula of Ray and Singer and the above computation yields: 

Theorem 10.3. For each finite CW complex, and each choice of skeleton, say X( d \ 


[Ut\o 
= nio 


(- 1 )' 12 

X( d ),i -I 

[det(d i+ 1 


• d*\image of d i+1 ) ■ (H x ,i) 2 ] { 1)! • 


This theorem gives a sequence of combinatorial identities since the terms det(d i+1 -d*{image of d i+ i) 
are expressed by Lyons’s formulae combinatorially. 

The papers of M. Catanzaro, V. Chernyak, and J. Klein, 013 E] also study aspects of the 
Reidemeister-Fraz torion in the combinatorial context. Their results are complementary to 
those given here. 
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11. Variant Settings and Examples 

11.1. Weights. One may easily adapt the above discussion to the context of cell complexes 
with weighted cells. Let {e,^} be the i-dimensional cells of a cell complex X and kR : 

R + be a positive real valued function. This weighting defines using the cellular 
basis an isomorphism, W l : Ci(X]R) —* Ci(X;R) multiplying a cell by its weight. Some 
familiar examples for graphs are: 

• One often introduces multiplicities for the edges. 

• For the normalized reduced graph Laplacian a vertex v is weighted by one divided by 
the degree of the vertex, see F. R. K. Chung HU- 

• In the classical Kirchoff circuit setting, the resistances of the wire edges are used as 
the weights. 

The paper of M. Catanzaro, V. Chernyak, and J. Klein [6j generalizes classical results of 
Kirchhoff on spanning trees, see also [4], to circuits, voltages, resistances, and is written in 
a weighted higher dimensional setting. 

Now all the present results readily extend to the weighted context. There are two formally 
equivalent methods to carry this out. 

In the first method: keeping the the standard cellular inner product, one replaces the map¬ 
ping d : Ci + i(X] R ) —* Ci(X; R ) represented by the matrix A and its adjoint d* : Ci(X ; R) —* 
C i+ i(X; R) represented by the matrix A by the mappings Cw = (khj) _1 / 2 <9(kkj + i) 1//2 repre¬ 
sented by (IR) -1 / 2 • A ■ (IFj+i) 1 / 2 and its adjoint £* = (kh^ • d* ■ (khj) _1//2 represented 
by (wHi ■ A ■ (IR) -1 / 2 . That is, one replaces the reduced Laplacian d ■ d* represented 
by A ■ A 1 by £ ■ £* represented by ((kR)" 1/2 • A ■ (W i+l ) 1 / 2 ) ■ {{W i+l f/ 2 ■ A* • (IR)- 1/2 ) = 
{Wi )- 1/2 ■ A-W i+1 - A t ■ {Wi)- 1 / 2 . 

In this approach the operators are changed, but the inner product is left alone. All 
covolumes are computed for the new operators with the old inner product. With these slight 
changes, all theorems can be recast in this extended context. 

In the second alternative approach: One recognizes that the characteristic polynomial 
of (Wi)~ 1 / 2 ■ A ■ W i+ 1 • A* ■ (IFj) -1 / 2 equals that of A ■ W l+1 ■ A t ■ (kR) -1 . Now if takes 
the adjoint of d represented by A with respect to the modified inner product < .. >w 
given by < e^e^fc >w— $j,k kF(ejj) then the adjoint, called here 9*’" is represented by 
Wi + 1 • A t ■ (kkj) -1 . That is, one may rephrase the above theorem in terms of d ■ by 
computing the covolumes utilizing the new inner product < .,>w but with d unchanged. 

All the results of the present paper can be readily extended to these weighted settings. 
For example, the analogue of Theorem 12.11 is 

Theorem 11.1. Let X be a d-dimensional cell complex. 

The k th elementary symmetric polynomial in the eigenvalues of the reduced weighted Lapla¬ 
cian C ■ £* are combinatorially expressed in term of covolumes of lattices with respect to the 
weighted metric < ., »w as follows: 

(Tfc(A) = [covolume{lattice image d* R (Cd{X] Z)) C Cd+i(X; Z) in Cd+i{X\ R))] 2 

^forests F of X of size k {^-‘coforests K of size k in forest F 

[i covolume {lattice image Bd-i(F] Z ) —> Cd-i{K ; R))] 2 ) 
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where the mapping B d _\(F\ Z) —> Cd-i(K; R) proceeds by inclusion into Cd-i(F] R) and 
then restriction to C d -i(K ; R) 

Here if there are no forests, or coforests of the required type, then that term is to be 
interpreted as zero. 

And the analogue of the corollary of the Trent type theorem in the weighted setting is: 
Let all covolumes in the next statement be computed with respect to the weighted inner 
product < . >w- Then an inductive calculation of the covolumes of Z^X; Z) C C;(A"; R) 

computed with respect to the weighted inner product is: 

Corollary 11.2. 

[covolume of lattice Z d (X]Z) C Cd(X; R))} 2 
= [covolume of lattice Z d -i(X; Z) C C d -i(X; R))] 2 ■ 

Wk,x ) 2 • [cavolume(B k (X;Z) C C k (X ; R))} 2 • (S TeTfc+1 (|^) 2 ) . 

We leave as an exercise to the reader the rewriting of the other theorems above into this 
enlarged setting. 

Remark: Likewise, all the results of this paper can also be rewritten to encompass weights 
which are allowed to be non-zero real numbers (or even non-zero complex numbers), e.g., 
using the normal principal values for square roots. An interesting example of positive and 
negative edge weights ±1 occurs in the work of K. Marasugi [23] on polynomial invariants 
of knots. 

11.2. Twisted Coefficients. Historically, Reidemciser-Franz torsion was introduced as an 
invariant of finite cell complexes with [/( 1) representations of the fundamental group, equiv¬ 
alently flat [/(1) bundes. The definitions were later extended to arbitrary unitary represen¬ 
tations. 

From this point of view the above combinatorial results are in the setting of the trivial 
U( 1) representation. As also done for some results on spanning trees by M. Catanzaro, V. 
Chernyak, and J. Klein [5j, in this section the above combinatorial results are extended to 
U( 1) representations. 

Given a representation p : 7Ti(X) —> U(l) of the fundamental group of a connected cell 
complex X into U( 1), one may form the chains with values in C twisted by the flat repre¬ 
sentation p, C*(X; E p ) together with it boundary mapping 

d p :C i+1 (X-E p )^C i (X,E p ) . 

Since p is unitary, the standard Hcrmitian pairing < z,w >K > zw on C extends to a Her- 
mitian inner product on i-chain with values in the local system E p : 

< ., >: Ci(x ; E p ) x C',;(Aq E p ) — > C 
Hence, the adjoint d *; Ci{X\ E p ) —» C l+x {X', E p ) is well defined. 

Now suppose that X is a d-dimensional cell complex with i-dimensional cells {e^}. One 
may consider the one dimensional vector space Ci(X\ E p ) of chain with support in that cell, 
say 

V[e.ij\ = {c G Ci(X] E p ) | c vanishes on cells not equal to e t j} . 
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This provides a canonical orthogonal direct sum decomposition of the i — chains : 

c,(A'l E„) = ® v[e,j] 

j 

For example, for any subcomplex S <Z X containing say k cells of dimension d and all of 
the d — 1 skeleton, one may consider the vector space of dimension k of C d (X; E p ) define by 

[S'] = {c G Ci(X; E p ) | c vanishes on cells not equal to e^} = V[ e d,j\ ■ 

ed,j cl d—cell of S 

Naturally there are boundary mapping (d\S) : —> V^_i[S']. 

Now picking bases elements of length one for the one dimensional vector spaces V[e d j] 
and V[e d -i,k\ the boundary mapping 

d:C d (X-E p )^C d ^(X-E p ) 

is represented by a matrix, say A, with rows and columns indexed by the d — 1-cells of A" and 
the d— 1 cells of A" respectively, or alternatively the one dimensional blocks, V[e d - gfc], V[e< 2 j]], 
respectively. 

Consequently, the algebraic lemmas 13.1113.21 apply, giving by the usual method 


a k (d p • d; : Cd-i{X] E p ) —>■ C d -i{ X; E p )) 

= a k (d; ■ dp : C d (x- Ep) -> C d (x- E p )) 

^ Blocks of size k y[ed j i 1 ]0-”®V[ed J i fc ] det(d* ■ dp\(V[e d)h ] ® ■ ■ ■ ® V[e d , ik ]) . 


Analogous to the above, dehne a subcomplex F C X to be a d-dimensional forest of size 
k for p if F includes the d — 1 skeleton of X, F contains exactly k cells of dimension d and 
the boundary mapping 

(d\S):C d (S ] Ep)^C d -i(S;Ep) 

is one to one. 

By application of the algebraic identity above, one gets the analogue of theorem 12.11 in 
this twisted context: 


a k (dp ■ d; : Cd-fX- Ep) —>■ Cd-i(X- E p )) 

^-‘d—dimensional forests F of size k det(Ap • Ap) 


d\S 

where Ap is the composite C d (F ; E p ) -b C d ~i(S; E p ) C C d - i(A; E p ). 

Of course, there are analogous theorems about forests and coforests and their dual notions. 
As mentioned above, the spanning tree case of coefficients in a fiat complex line bundle 
where k = dim B d _fX] E p ), is discussed in M. Catanzaro, V. Chernyak, and J. Klein [5] 
and there related directly to circuit theory. 


11.3. Manifolds. When K is a triangulation of a closed n-manifold, one may introduce the 
dual cell subdivision, K * whose (n — i)-cells, D(e it j), are indexed by the baricenters, say 
<7 ij of the i-dimensional simplices a hr The dual cell D(aij) consists of the sets of simplices 
< & % n , (Tjp i j 2 , • ■ • ,a n j n _ i+1 > of the baricentric subdivision of K where the baricenters 
correspond to an increasing flag of simplicies of K starting with a,j = a hn . That is, the 
simplicies a i+k -ij k of K are of dimension i — 1 + k, with k = 1, • • • , n — i + 1 and a % ,j = 
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a i j 1 C cri+i,j 2 C ■ • • C cr n j n _. +1 . The bijection <Jij t->- establishes an isomorphism of 

chains 

Ci(K; Z ) = C n _i(K] Z) and C^K; R ) = C n -i(K; R ) 

under which the boundary d of K and its dual d the “coboundary” for the cellular in¬ 
ner product on K , are identified with the coboundary for K * and the boundary for K *, 
respectively. 

Hence, in this special context, forests of K correspond to dual forests of K *, coforests of 
K correspond to dual coforests of K *, etc. 

12. Enumerative combinatorics of the simplex following Kalai p3] 

Let A n _! denote the standard n — 1 simplex with n vertices, say x±,--- ,x n in this order. 
Regard A n _! as the join 

A n _i X\ k A n _2 

where A n _ 2 is the standard n — 2 simplex with vertices x 2 , • • • ,x n . For a simplex, say s in 
A n _i, let s have its inherited orientation from this choice of ordering; for emphasis, denote 
the oriented simplex s. As s ranges over the j-simplices of A n _i, respectively A n _ 2 , the 
elements 1 • s provide an integral basis for the j-chains, Cj(A n _i, Z), Cj(A n _ 2 ; Z). With 
respect to these cellular choices of basis, the boundary maps 

■ Ck{ A„_i; Z) —>■ Cfc_i(A n _i; Z ) and <9a„_ 2 : Ck{ A n _ 2 ; Z) —> C y fc_i(A ri _ 2 ; Z) 

are represented by the respective incidence matrices, and have adjoints, with respect to the 
cellular bases, d\ n l and represented by the transposes of these incidence matrices 

respectively. 


Note that the integral d-chains, Cd( A n _p Z), for d < n — 1 have a spanning tree consisting 
of all d-cells, say s, containing the initial vertex X\. Moreover, associated to this spanning tree 
the integral d-cycles, Z^(A n _i;Z) for d < n — 1 have a basis consisting of the boundaries 
d(xi -k t) where t ranges over the d-cells in A n _ 2 , that is the d-cells not containing x\. 
Consequently, the above theorems apply in this context. 


Now consider the standard cellular inclusion : 

i '■ Ck~ i(A n _ 2 ; Z) — y Ck- i(A ri _ 1 ; Z) 

In the cellular basis, the adjoint is just the cellular restriction mapping 

n (A V\ restriction = i* n / * 

V'fc- l\LA n -li^) —r Ofc_i( k ZA rl _ 2 , Z J . 

Following, Kalai, consider the composite : 

: Ck( A„_i; Z) — > Ck- i(A n _ 2 ; Z) 

It is represented by the restricted incidence matrix, say /*, obtained by deleting the rows 
corresponding to k — 1 simplices not in A n _ 2 . In these terms, one of Kalai’s results takes 
the following form: 
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Theorem 12.1 (Kalai ). The composite 

{i* • <9a„_i) • (** • ^a„_i)* : Ck- i(A n _ 2 ; Z) —> C k ~ i(A n _ 2 ; Z) , 

which is represented by I* ■ (Iff, has eigenvalues 1 and n with respective multiplicities (”“ 
and ( k )• 

Note that, as expected, the sum of these multiplicities is ( n “ 1 ) = rank Cfc_i(A n _ 2 ; Z). 


D 


The result easily gives the eigenvalues and their multiplicities of the Kirchhoff matrix 
representing the combinatorial Laplacian 


<9 An _ 2 ' ^a„_ 2 • Cfc- i(A„_2 : Z ) —> C fc -i(A n _2 : Z) 

as follows: 

For s a k — 1 simplex in A n _ 2 , one computes : 


(** ' d An _J • (i* • 9 An _ 1 )*(l ■ s) = (i* ■ d An _J • dA n _! •*(!■«) = 
= (** • 9a„_J ■ (1 -a:i*3 + i ■ (<9 An _ 2 (l ' «))) = i*( 1 ■ s-Zi * 
= l-5 + a An _ 2 (^ n _ 2 (l-S ) 


d An _ 2 (!-s ) + dA n _ 2 (dl n _ 2 


( 1 - 5 ))) 


That is, one sees that 


(** ' dA^) • (i* • dA^)* = Id + d An _ 2 (dl n J 

yielding the following corollary to Kalai’s result. Let Bd k be the matrix representing the 
boundary mapping <9 A „_ 2 : C k { A n _ 2 : Z) ->■ C fc -i(A n _ 2 : A). 


Corollary 12.2. The combinatorial Laplacian 

d An _ 2 • ^ A „_ 2 : C'fe-i(A n -2 : Z) —» Cfc_i(A n _ 2 : Z) 

is represented by the matrix Bd k ■ ( Bdf) 1 and has eigenvalues 0 and n — 1 with respective 
multiplicities (£“ 2 ) and (”“ 2 ). 


Now there is a natural integral and geometric basis for the integral k — 1 cycles of A n _i, 
Zfc_ 1 (A n _ 1 ; Z), provided by the following map from k — 1 chains on A„_ 2 to k — 1 cycles 
which is an isomorphism: 

^ : C'fc-i(A n _ 2 ; Z) —> Z k -i( A n _!; Z) 

1 • s (->• <9 An _ 1 (xi *s) for s a k simplex in A n _ 2 

Since S An _ 1 (xi*s) = 1 • s — xi*d An _ 2 (l • s), the composite i*<h(l-s) equal just 1-s. That 
is, i* ■ <f> = Identity. Hence, <f> is one to one. 

Next given a k — 1 cycle, z G Z k -i(A n ^ 1 ; Z), then z = <3>(i*(A)) + (z — <3>(i*(A)) where 
i*(z — $(i*(z)) = 0. Hence, (z — <f>(i*(A)) = X\ -kc for some k — 2 chain c supported in A„_ 2 . 
But then 0 = i*(<9 An _ 1 (A — < h(i*(^))) = i*( 9 An _ 1 (xi * c) ) = i*( c - aq * <9 An _ 2 (c) ) = c. This 
proves that c = 0 and so that z = $(i*(z)) when z is an integral k — 1 cycle. 

In toto, this shows that is an isomorphism. 

In particular, the cellular basis for C k ~ i(A n _ 2 ,; Z) gives the geometric basis {$(1 • •?)} for 
the k — 1 cycles. As expected, using the Kalai result again, one may use this method to find 
the eigenvalues and multiplicities of the mesh matrix for k — 1 cycles. 
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Let J k 1 represent the inclusion Z fe _ 1 (A n _ 1 ;Z) C C'fc_i(A n _ 1 ; Z) with respect to the 
geometrical basis {$( 5 )} of k — 1 cycles and the cellular basis of k — 1 chains. That is, J k ~ l 

represents the composite C'fc_i(A n _ 2 ; Z) A Z^~ i(A„_i; Z) -b Cfc_i(A n _i; Z) with respect to 
the cellular bases. 

Corollary 12.3. For the integral basis of Zk-i(A n ~i, Z) provided by the geometric integral 
basis <h(l • s) where s ranges over the k — 1 simplices of A„_ 2 , the mesh matrix for cycles 
representing the pairing 

Zk-i(A n -u Z) x Z/ : -i(A n _ 1 ; Z) —> Z 

with respect to this basis, or equivalently of the matrix has eigenvalues only 

equal to n and 1 with respective multiplicities (”“ 2 ) and ( n , 2 ). 

Since Z^(A n : Z) = Bd( A n ; Z), the cycle and boundary mesh matrices are identical in this 
case. 

To prove this it suffices to compute and simplify the entries with s, s* two k — 1 simplicies 
in A n _2 : 

< <h(l • s), <h(l • s* >=< d A (x\ -k s), d A (xi * s#) > 

=< 1 • s-si*3 a „_ 2 (s1 , 1 • s* - x x k d An _ 2 (s*) >= S S:S #+ < d An _ 2 (s),d An _ 2 (F&) > 

— < s,9 a„_ 2 d An _ 2 {s^) > 

Now by Lemma 1X21 the non-zero eigenvalues of d An _ 2 d An _ 2 acting on Ck-i(A n _ 2 ; Z) are 
the same as that of the combinatorial Laplacian, or Kirchhoff map, d An _ 2 d An 2 acting on 
Cfc_ 2 (A n _ 2 ; Z). By corollary 112.21 this last has a single non vanishing eigenvalue, narneby 
n — l, of multiplicity (”“ 2 ). Hence, the mesh matrix above has an eigenvalue n of multiplicity 
(£ Zi ) an d eigenvalue 1 of multiplicity ("^i 1 ) — {^Zl) — ("fc 2 ) as d a i me d. 


Kalai has given an additional refinement of his result recorded here in theorem 112.11 
For each vertex Xj introduce a real positive indeterminant, say aj > 0. Introduce a weight 
for each face / =< ay 2 , • ■ • , x H > with w(f) equal to the product a tl • a, 2 • • • a*, of the 
corresponding indeterminants of the vertices of the face /. This defines an inner product on 
each of the j chains for which the elements {1 • // \Au(/)} are an orthonormal basis. These 
are called the normalized cellular basis, Let D(k,m ) be the diagonal self mapping of the k 
chains on A m which multiplies a k face / by w(f). For convenience, also let D(k,m ) also 
denote the diagonal matrix which represents D(k,m) in the cellular basis. 

Now the Kalai mapping i* ■ d An _ 1 : Cfc(A n _i; Z) —* C'fc_i(A n _ 2 ; Z) with respect to this 
new normalized cellular basis gives a reduced incidence matrix, say (/ r fc )^ which is related 
to the original reduced incidence matrix if by 

(If)# = D(k-l,n- 2)" 1/2 If ■ D(k , n - 1) 1/2 . 

Kalai’s extension states that the matrix 

(If)* • ((. If)*y = D(k -1,7 1 - 2) -1 / 2 If -D(k,n-1) (ifY ■ D(k , n - 2 )~ 1/2 
has eigenvalues cq and E" =1 Oj with multiplicities (ffZi) an d (”^ 2 ) P3]- 
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Correspondingly, the boundary mapping <9 a„_ 2 : Ck- i(A n _ 2 : Z) —> Ck- 2 (A„_ 2 : Z) with 
respect to the new normalized cellular basis is a matrix, say [Bd k ~ 1 )^. It is related to the 
standard boundary mapping Bd k ~ l by 

(. Bd fc “ 1 ) # = D(k 2) -1/2 (5d fc_1 ) £>(fc, n - 1) 1/2 . 

The normalized Kirchhoff mapping is then represented by 

(. Bd k ~ l )* (( Bd = D(k-l,n-2)~ 1/2 (Bd k ~ l ) D(k,n-l){Bd k ~ 1 ) t D(k-l,n-2)~ 1/2 . 

A slight extension of the above argument shows that this matrix has eigenvalues 0 and 
E cij with multiplicities with respective multiplicities (”“ 2 ) and ( n “ 2 ). This follows from 
Kalai’s results by the method used in the unnormalized case above. 


Correspondingly, under the identification $ : Cfc-i(A„_ 2 ; Z) —> Zfc_i(A n _i; Z), one has a 
bilinear “mesh” pairing 

Cfc-i(A n _ 2 ; Z) x C'fc_i(A re _ 2 ; Z) —> Zk- i(A n _i; Z) x Z*_i(A n _i; Z —> R 

where the last bilinear pairing used the chosen weights w(f). 

For the choice of normalized cellular basis for Cfc_i(A„_ 2 ; Z), it is easily seen, by the 
method above, that this matrix has eigenvalues (i\ and E^R 1 a 3 of respective multiplicities 
“’) and ( n “ 2 ). This follows from Kalai’s results by the method used in the unnormalized 
case above. 

Extending in several directions some of the formulae of Kalai, the complete calculations in 
this section of the eigenvalues of the Kirchhoff matrix and of the two mesh matrices in each 
dimension of simplices, yield, by results developed in the preceding sections, many combina¬ 
torial identities involving appropriate subobjects of simplices for each of the coefficients of 
the corresponding characteristic polynomials. 
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